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Abstract 

The serial harness introduced by Hammersley [7] is equivalent, in the 
Gaussian case, to the Gaussian Solid-On-Solid interface model with par- 
allel heat bath dynamics. Here we consider sub-lattice parallel dynamics, 
and give exact results about relaxation dynamics, based on the equiva- 
lence to the infinite time limit of a time periodic random field. We also 
give a numerical comparison to the harness process in continuous time 
studied by Hsiao [8] and by Ferrari, Niederhauser and Pechersky [4j[5]. 
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1 Introduction 

Let L be a positive even integer and let the initial condition hP = {h® : i <G 
Z/LZ} be distributed according to the un- normalized measure 

Md/ i )=IIe"" (/l?+1_ ^ )2 d^ (1) 

i 

where dft.^ is the Lebesgue measure over R. The index i runs over Z/LZ, which 
corresponds to periodic boundary conditions. The measure (JTJ) may be consid- 
ered as a finite volume Gibbs measure with Hamiltonian 



ff(i. )4E(t-ft.°) 



2 

t 

The corresponding sub-lattice parallel heat bath dynamics is defined by 

■(hi-H h i-i + h i+l)) 2 



dh t \h t - i )= n 

in 

• LI m-hl^Hdhl/noim. (2) 



c 

i-\-t even 



i-\-t odd 



l 



where the normalization of the probability is a finite constant, independent 
of ft, 4-1 . The stochastic process defined by ([2]) is intermediate between Ham- 
mersley's original serial harness [7] and the harness process in continuous time 
[SI SJ [5] . Various sub-lattice parallel stochastic dynamics for interface models 
have been studied, e.g. in [TJ [3J, showing a closer similarity with continuous 
time dynamics than with fully parallel dynamics. 

The heat bath dynamics leaves invariant the Gibbs measure which motivated 

it: 

' n{dh t - 1 )P( dh* I h 1 - 1 ) = M (d/i 4 ) (3) 

As the initial condition h is already distributed with the measure \x, we have 
a stationary problem. Our main result is a computation of space-time corre- 
lations, in the thermodynamic limit L — > oo. The correlation function of two 
space gradients at time and space separation (2t,j) will be denoted gn(t,j), the 
correlation function of two time gradients at time and space separation (2t,j) 
will be denoted g 22 (t,j). The time separation 2t corresponds to t updates at 
each site between the two events: 



g n (t,j) = Km E(hf +2 - hf){h° 2 - h° ) , t>0 (4) 

g 22 (t,j) = lim E(hf +2 - hf)(hl -hi), t>0 (5) 

and similarly 

g 12 (t,j) = lim E(hf +1 - hf_M ~ &o) . * > 1 (6) 

g 21 (t,j) = lim E(ftg* - h%- 2 ){h? j+l - h ^) , t > 1 (7) 

Proposition 1. Let /ii 0,t l be distributed according to {IPfiP, /or eac/i t € Z+. 
T/ien /or eac/i t G Z + and j e Z i/ie limits (0) (EP ® (Q' e:E *si arcd satisfy 

gu(t,j) = ; if < j < 2t , (8) 

V ' (t -i)\{t + i)\ -■> - i \ > 

gu(t,-j) = -gii(t,j) if j>l (9) 

= if |j|>2*>2. (10) 

ff22(t,i) = -i[sn(t-i,UI)-flii(t,|i|)] if t>\ (11) 

gi2(t,j) = -gi2{t,-j) = -g 21 (t,j)=g 21 {t,-j) (12) 

= 5n(i-l,|j + l|)-9ii(t-l,|i-l|) (13) 

Moreover, as t —¥ 00, uniformly in j G Z, 

flll (t,j) = J=e"it +0(t- 2 ) (14) 
v 7ri 
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922 (t,j) 



1 



4Wirt 



(l-£)e-£+0(t-') 



lim E(h) - /i-)(/iq 



9l2(t,j) 



V 7T 



2j _3l 
e 4t 



2t 



Vt 



— 7= I due 



-u 2 /2 



(15) 
(16) 

■OQnt) (17) 



Remark 1. [9]) l 1 (J\ ) fLlfy follow respectively from space symmetry, causality, and 
the detailed balance condition. More identities can be extracted from (0) and also 
from the loop condition ( the sum of gradients around a closed loop is identically 
zero). In particular gii(f, 0) = kg\%{t, 1) Vt. 

Remark 2. Proposition^ conveys information for \j\ <C Vthit. 

Proposition Q] is proven in Section [5j It is based on the equivalence in law of 
the space-time field fti '*' with the infinite time limit of a space and time periodic 
random field, which is naturally diagonalized by Fourier transform. This random 
field is defined in Section [2j The Fourier transform diagonalization is performed 
in Section G3 The proof of equivalence is completed in Section |H Generalization 
to arbitrary dimension is outlined in Section [6] A numerical comparison to the 
harness process in continuous time is given in Section [7] 



2 Space-time periodic field 

For T a positive even integer, the marginal space time field 

h={h\ : (i, i) € ( {0, . . .T - 1} X Z/LZ ) n {t + i even} j (18) 
is easily checked to be distributed according to the un-normalized measure 

TV, s / TT -(^-KCl+Cl)) 2 , 

niuah)={ n e dh * 

t-\-i even 

[] e"^"^ 2 ^ ) (19) 



where "free" refers to the time T — 1 final condition and the range of t in the 
product is 1 < t < T — 1. The corresponding "space-time Hamiltonian" is 



*s = E ( h *~ \ h t\ - l'r + i) 2 +\ E ( h °- h Q 2 ( 2 °) 

i+t even i CVC11 

1<£<T — 1 



A good feature of nj r ^ c is that its marginal at time t is known exactly. However, 
in order to compute time correlations by Fourier transform, we are going to use 
periodic boundary conditions in the time variable also: let 

h = | h\ ; (t,i) e (Z/TZ x Z/LZ) n {t + i even} j (21) 
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be distributed according to the un-normalized measure 

&h)= II e dh > ( 22 ) 

i+i even 

The corresponding "space-time Hamiltonian" is 

1 1 \ 2 

■i+i even 



z+t even 

E fa - ^=r ~ \^ + \) 2 + E fa ~ ^f-i 1 - ^ 



■i+teven j even 

KKT-1 



(23) 




Fig. 1: Even space-time sub-lati 



The last term in (|23]l or ([20| is necessary in order to have only one mode 
distributed according to the Lebesgue measure, for uniform global translations 
of the system. Fig. 1 shows in solid line the interaction terms common to (|23[) 
and (|20l) . and, in dashed line, the interaction terms corresponding to the last 
term in ([2"3"| or (|2Hj) . 

Proposition 2. Let L,T\,T be positive even integers. Let h^^. and hj^ be 
random fields distributed according to fi^^. and (J-J^ respectively. Then, as 
T — > oo, the marginal ^per|o<t<T -l conver 9 es * n distribution to hj^. In par- 
ticular, the one-time marginal h~, er converges in distribution to hf Tee , distributed 



4 



according to fi(dh°). And, extending the random fields to the full lattice with 
(hper)j = (^pcr)i f or t + 3 °dd, and similarly for h( Tee , then h® el converges in 
distribution to hf Tee , distributed according to fi(dh ). 

Proof. Each random field has one real component distributed according to the 
Lebesgue measure, the same for all random fields. We need only consider the 
gradient fields. The statements about the one-time marginals follow from the 
convergence of the one-time covariance matrix, e.g. E(/j? — — hj_^) 

for i,j — 1, . . . ,L — 1, which are linear combinations of E(/i° — ft.§) 2 for j 
\ ..... I. 1. This is a computation, given in Section HI We thus have, for the 
gradient fields, 

P^(dh°) — >P&ee(dfc ) as T^oo (24) 
On the other hand, 

m TLf dh [l, Tl -l] | ^0) = pi ee ( d/l [l.r 1 -l] | ^ h T = h 0j 



- per V 



5 £ee(d7i [1 ' Tl - 1] |/i°) as T^oo (25) 



Then, for the gradient fields, 



- free 



d/.! ^- 1 ] ) (26) 



□ 



3 Fourier transform 

In order to compute the Fourier transform of the space-time periodic Gaussian 
field, it is convenient to set 

h 2 2 l +1 =0 and h%&l = Q Vt,t (27) 

Then, for (i/, k) € (Z/TZ x Z/ZZ), 

k ~ ^fW 2^1 2-^ e n 3 ' ""i-fc ~ ' n k+L/2 ~ n k 

V L1 j =0 t =0 

(28) 

and 

1 -2m^--2m^ . 

v ljl k=Q u=0 
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and 



L-l T-l 



n pcr - E E ~ 2 ft *'-i ~~ 2^+^ 



j=0 t=0 

TT ^ 



/i£/i£,e 



-2i7r(fc + fc')i -2i7r(i/ + ^')T 



('- 
('- 



1 2i7r^ + 2i7rf i -2i7r^ + 2i7rf s 

— -e 

2 2/ 



1 2i7Tf +2i7T^ 1 ~2i7T% +2i7T^x 

7:^ - -e 

2 2/ 

1 2i7T^+2i7rf 1 -2i7T^ +2i7T^x 

7:6 - -e 

2 2/ 



/ 1 -2i74-2i7r£ l 2i7rf-2i7r^ x 



1 — 2 cos27T— cos27T— + cos 2ir 



i) 



(30) 



The Fourier transform may be cast into an orthogonal transformation of 
the LT/2 random variables h\ . Setting h k = a k + ib k , the new LT/2 random 
variables a k 's and b k 's are chosen as follows: the orbit of any (v, k) under 
possible combinations of 



0, k) -> (T - v, L - k) and (v, k) -> [y + T/2, k + L/2) 
has 4 elements, except for the 2-clement orbits 

{(0,0), (T/2, L/2)}, {(0, L/2), (T/2, 0)} 
{(T/4, T/4), (3T/4, 3T/4)} , {(T/4, 3T/4), (3T/4, T/4)} 
corresponding to real h k } s. Choosing one element per orbit, we get 



(31) 

(32) 
(33) 



>LT 



[a° + (-yd 



T/4 
3L/4 



L/2-1 



+ 



+ 



'LT 
4 



L ° fc L 

k=l 
3L/4-1 

E H 27r /7 + T)^ /4 + sin ( 27r T + T)^ 

iji fe=L/4+l 



E [cos27r^a£ + sin27r^6£ 

kj Tit, T/4 



- k 3 _l nt \ uT/i 



(34) 



T/4-1 L-l 



+ 7^f E E[ cos ( 27r Y + 27r ^^ + sin ^T + 2 ^f^^ 



i/=l fc=0 
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The jacobian of the transformation from h\ 's to a k 's and b k 's is actually 2 TZ 7 4 . 
The new measure is 

ex P (-^ e L r )n d <n d5 fe/ norm - ( 35 ) 

v.k v,k 

where the index sets for the a k 's and b k 's are as in the formula (fM|) for hj , with 
a total of LT/2, and 

H™ = 2 7^K) 2 +4 7^[K) 2 + (&D 2 ] (36) 

2— orbits 4— orbits 

with Y k defined in (00]), so that E|/i£| 2 = l/(4 7 £) for all (v, k) ^ (0,0) or 
(T/2,L/2). We have one zero mode a® — a^y 2 , distributed with Lebesgue 
measure, and soft modes, Gaussians of large variance, around the zero mode: 

Lemma 3. h® — h^j^ * s distributed according to the Lebesgue measure. All 
other h k are independent centred real or complex Gaussian variables with 

^\K? = - " — i r (37) 

4 1 - 2 cos 271-f; cos 2tt£ + cos 2 2tt| 

Eh k h k , is non zero only if (y, k) and (y 1 k') belong to the same orbit, with h k and 
h k , complex conjugate: (y 1 ', k') = (T — i>,L — k) or (y' ', k') = (T/2 — v, L/2 — k). 

4 Equal time covariance 

The equilibrium measure ([T]) can also be diagonalized by Fourier transform, 
which yields 

^-^^^ rs (38) 

u k ^ Q i cos L 

where k 6 Z/LZ. A change of summation index k — > h — k leads to equivalent 
formulas, differing slightly according to the parity of j. Averaging the two 
formulas yields 

1 — cos — - 



E(hS - hlf =i Y C °\, L k if J even 

y J IV l Z^f i _ 2 2£fe 



1 1 1-COS^COS^ 

- + - > n £ if J odd 



(39) 



Here of course L = oo would be simpler, with M(hj — ft.[j) 2 = j, but our aim is 
to complete the proof of Proposition 2, where L is finite. Let us now compute 
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the analogue for the space and time periodic field: 

1 _ 2l7T^ -2i7r^ 

E per (/i° - h° f = — Yl ( e i ~ 1 )( e L -1)EW (40) 
Therefore, with j even, using Lemma [31 

9 -2i7T-^ 2i7T^. 

E per (/^°) 2 = A^ (c * _i)( e * -1)EK| 2 

1 1-cos2tt^ 
"iT,,^,^ 1 - 2 cos 2tt^ cos 2tt£ + cos 2 2tt£ 



1 1 Z" 2 ^ l-cos2^M 

— > — duj r — r as T -5> co 

27r7 1 - 2coswcos27r# + cos 2 2tt£ 



L 1 wo *" L 

1 ^ l-co^irff 

i . ~ 1 - cos 2 2tt£ V 7 



fc^0,L/2 i 

where the last step comes from the identity [6], p 366 
1 f 27r , cosncj 



— / - = - (42) 

2ir Jr. 1 — 2ocoso; + a 2 1 — a 2 



The covariance (|4"Tj) is indeed the same as the covariance (f3"9")) . Similarly, for j 
odd, 



2'nrkj _ 2\irv 2'i-Kkj , 2i7rt/ 



2i7rfcj _ 2i7rt/ 2i7rfcj . 2i7rt/ 

1 _ 1 - cos(2ttM + 27rf ) 

IT ^ , „ 1 - 2 cos 2tt£ cos 2tt£ + cos 2 2vrf 

1T-1 1 ^ ^ 1 - cos(2tt^ + 27rf ) 

T T + LT ^ r ^1-2 cos 27rf cos 2tt£ + cos 2 2tt£ 



1 1 1 f 2 * 1 - cos(2ttM + w ) 

- + - > — / dw £ r as T ->■ oo 

I L u _£-f /o 27T / 1 - 2coscjcos27rf + cos 2 27rf 

-M 



MO.L/2— " u * £ '"i 



1 1 ^ 1 f 27T 1 -cos(2tt^-)cosw 



T T , ,„ 27r Jo 1 - 2coswcos27r-| + cos 2 27r-| 

1 1 v-> 1 -cos27r^cos27r| 
= I + I fe it/2 l-cos 2 2^ 

(43) 

which again is the same as (f39|). This completes the proof of Proposition 2. 



5 Proof of Proposition [T] 

Proof of 14)); using ([29)) and Lemma [3) with t even, 
E per (/i* +2 -/i*.)(^ - fcg) = 



2 v.— — -x-^..^ 



^E( e L T - e L T )( e L -«i 2 



LT^ 1 



cos 2ir^ cos 27rf (l - cos 4tt£) 



£T J^»l- 2 cos 2tt£ cos 2tt# + cos 2 2tt£ 

f 27T ^ cos 4>j cos wt (1 — cos 2c 
47r 2 J "~ 7 1 — 2 cos u; cos + cos 2 



27r j I -2 * j- cos^'coswi(l-cos20) 
dw / d<f>- — ^— £j asi,T^oo (44) 
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so that 

1 /" 27r j /"^j cos <j>j cos ut (1 - cos 20) 



lim E(/iL 2 - - = J- / dw / ., , 

l^oo v JM ^ uy 4?r 2 7 7 l-2coswcos< 

1 

7T 

2- t+1 i! 



2tt 

t 



I) 



cos </>j (cos 0) 



2 • 2 ■ 

(45) 

where we used (|42|) from the first to the second line, and the last line assumes 
j < t, otherwise the result is zero. Changing t into 2t yields (|SJ). Stirling's 
formula then leads to (11411 . 



Proof of (ElPmP: 

E per (/ l *+ 2 -/i*)(^-/ l 0) = 

9 -2i7r^. - 2i7r t/(t+2) -2i7r^ - 2i7r— 4i-7r- 

= Z^zJ e " e )( e -W^fel 



— y y- 

LT ^ V 1 



cos 2tt^ cos 2tt^ (1 - cos 4tt^) 



, . 2cos2tt£ cos2tt£ + cos 2 2?rf 

I/jfcQ, T/2 K 1 L L 

1 /" 27r , /" 27r , , cos0jcoswt(l-cos2w) 
^ Jo Jo 1 — 2 COS COS (p + COS (/> 

so that 

1 /" 2,r , /" 27r cos 0j cos (1 — cos 2w) 



1 

47T 



1—2 cos w cos + cos 
cos^j (cos^) (t_2) - (cos^)' 



(47) 



where (|42p was used once more. Comparing with the second line of (|45[) gives 
(fill) , which combined with (fl"4"|) gives (ITS1) . 

Proof of FIB: 



E pei (h t j+1 -h t j _ 1 )(h 2 -h° ) = 

2 „ -2i7rf - 2i7rM 4i7r^ . . -2i7r4 -2i7r£ N A , 

4 x-^ — 2i7r^ 4i?r^ . . . , A 

= -^E e (e - 1 )«n2^sin2 7 r^E|^| 2 (48) 
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where the space symmetry j — > —j was used to get the last line. In order to 
use time reversal symmetry, coming with the detailed balance equation, let us 
compute similarly 



4 

LT 



— E c ( e -l)sin2^ sin2vr feE|^| 2 (49) 



Time reversal shows that f|49|) is the opposite of (1481) , which therefore equals 
E pei (h t j+1 -h t j _ 1 )(h 2 -h° ) = 
= j^J2( cos2n T -cos27r^^)sin27r^ sin27r^E|^| 2 

1 /" 27r /" 27r (coswi-cosw(t-2))sin0jsin0 
— > — r / dcj / d0 — ' 5— as L, T ->• oo 

7T ,/q Jq 1 — 2 COS CJ COS + COS^ <p 

(50) 

so that, using (|42]) once more, 

1 Z' 277 

lim E(/i* -/i*" 2 )(/in-/io) = - / d0(cos< ? !>)(*- 2 )(cos(/)(j + l)-cos0(j-l)) 

(51) 

Comparing with the second line of (|45|) gives (p~3|) . which combined with (| 14[) 
gives (HI 



Second proof of iu5j) : We give here a second proof of (|T5|) . not relying upon 
CD). Let 

g{ U ,j) = (1 - cos 2w ) J T - i ___^_-^ = g{—u,j) = - W> j) 

(52) 



so that 



/" 27r , /" 27r , , cos <f>j cos (1 - cos 2w) 



1 — 2 cos w cos 4> + cos 2 1 



2tt 



^ /-27T 



7T/2 



= -~ / d W e-V( W ,j)+cc. (53) 
lj Jo 

In order to perform the integral over <f> and estimate g'{ui,j) near cj = 0, we 
decompose 

1 1 1 



1 — 2 cos uj cos + cos 2 6 1 — e 2lw 1— e lu cos 



(54) 
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and use the residue theorem: 

i_ ff l-e lw cos0 u ' J_, 



(j) 2 -(2- 2e-^) 

= Otr 1 ) + 2e~ iu , 2l?r — e 1 ^ 2 ^' 
2V2 - 2e- i <" 

- ©(r 1 ) + 0(^ 1/2 ) + ^vre^- 1 / 2 e -v^e--/ 4 W 1 /2 J 
* d(j) cos $7 



1 — e lw cos < 



(55) 



because the integral with e -1 ™, with the contour closed in the lower complex 
half-plane, gives an equal contribution. The 0(j ) is a regular function of oj, 
bounded by const, j as j — > oo. Therefore 

j) = O^r 1 ) + 0{^' 2 ) + V^ne-'^W/ 2 G -^~' m/ ^ 1/2 3 + c . c . ( 56 ) 

and 

g>(u,,j) = Oir 1 ) + 0{^' 2 ) + iTrj e-^e-W^Vaj 

+ JLe--/ 4 ^- 1 / 2 e-^e-^/ 4 ^ 1 / 2 ^ + c c (57) 
V2 

In (|57p. the phase factor oscillating with cj is exp(ia; 1 / 2 j) in both terms shown 
in full, whereas the c.c. counterparts will have exp(— iui 1 / 2 ]). When integrat- 
ing against exp(iu;i) the leading contributions will therefore come from the c.c. 
counterparts, which will give a stationary phase region. This agrees with the 
result of the calculation below. 

7r ' 2 • * AToW4, ,1/2,; , n? ■ f°° ■ t A? iw/4. ,1/2 „■ 

ii A) 

-i7r/4 

/•e oo _ 

= 0(t- 1 e~V* i ) + 2i / xdxe~ tx2 - iV2:ix 
Jo 

(■OO 

= 0(t- 1 e"~vf J ') + 2i / xdxe- tx2 - iV2]x (58) 
Jo 



r/2 d WW - 1 / 2 c i - t -^e i ^ 1 / 2 ^ or i )+ r^-xfijut-yft&Wj 

Jo 



PC ' OO 

= 0{r 1 ) + 2e i7l/4 dxe^ 2 -^* 
Jo 

/>oo 

= 0(r 1 ) + 2e i7l/4 / da; e"** 2 - 1 ^ 
Jo 



(59) 
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so that 



/(t,j)=Or 2 )+i— [°° xdxe- tx *-^ x -^= [ dxe-^-^+c.c. 
* Jo W 2 Jo 

(60) 



Then 



* io * Lao %/2t 5 / 2 

(61) 

and 

= / dxe tx W2jx +c.c.= / die fI W2jx = = — -e at 

(62) 

Altogether 

47T 2 / 7 2 N f. 

f(t, j) = 0(t- 2 ) - =— (1 - ) e- St (63) 



V^f 3 / 2 V t 
which completes the second proof of (|T5|). 

Proof of {23). ■ 

n » n 2 -2i7rM_2i^ -2i7rM 

lEperC^j - h°)(hl-hl) = — £(e T - e ) • 



(e T -1)E|^| 2 



^ 1 - 2 cos 2tt£ cos 2tt| + cos 2 2tt| 

1 f 277 f 271 (l - COSUjt) COS 07 

— / dw / d0 i '- -Vt ^ T -> oo 

4tt 7o Jo 1 — 2cosw cos</> + cos^ </> 



1 ^ (cos2ttM)(i - cos27rf) 



1 r' 2 f 27T (1 - coswt) COS07 
= — / dw / d<£— i- 7 W 2 - *,jeven 64 

7r JO JO 1 — 2 COS W COS (f> + COS 



The integral over <p is done like in the proof of (|15p . Then, using 



1 — cos uit = oj I dssinws. 



o 



E{h] - h°)(h - h Q Q ) = 0(lnt)+ 

+ _L_ /" ds T /2 du e -W4 w -i/a sines e -^ e "" /4wl/2 J + c.c. (65) 
2 ' 7r j Jo 
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Then, setting lu — \x 2 and using Cauchy's theorem, 

/>e — ^ oo 



- dtdc-V^^-^ 4W5 J =0(«-i) + i / dxe" 
2i Jo 1 Jo 

= 0(s- x ) + - / dxe-™ 2 "^ 
1 Jo 



-V2jx 



(66) 



whose main part will cancel out with its complex conjugate. Similarly, setting 
u = —ix 2 and using Cauchy's theorem, 

/ " d W e-^^^e- iws -^ e l4w ^'=0( s - 1 )+ / dxe- s - 2+i ^ 
2i Jo Jo 



0(0 + I dxe- sx2+i ^ x 

(67) 



and 



3G 



+c.c. = / dxe-" aS+iV ^ a = \ -c--s (68) 

J-oo V s 



Altogether 

E(/>* - h^hl - h ° ) = -j=J*^ s e-fi + 0(ln t) 



7 e 2t u +0(lnt) 



7T ./i 

27 



e 



which completes the proof of Proposition [T] 



O(lnt) (69) 



6 Higher dimension 

In arbitrary d > 1, ft = {/i? : i G (Z/iZ) }, the initial measure 

-±(/i° - ^ 2 

M(d^= n c ik ( 7 °) 

it~ji=i * 

is invariant under the dynamics defined by 

-d( h l - y /i*-M 2 

p(dfc*ifc*- i )= n - J 



e 

even 



II - /x*- 1 ) IJ / norm. (71) 

|i|+t odd i ' 
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where runs over the 2c? neighbors of i. The space-time Hamiltonian on the 



d+ 1-dimensional torus (Z/TZ) x (Z/LZ) d is 
The Fourier transform is defined as in d = 1, using k = (fci . . . kd). Then 



(72) 



r i d k i d 



~L 



= El^| 2 7k 

i/,k 



kl 2 



(2ttz/) 2 f 1 (27rk) 2 ^ t 1 ^(27rk) 2 ^- 



T 2 



2d L 2 



d 2 



v 2L 2 ; 



(73) 



and the autocorrelation can be computed, yielding the expected ln(i) in d = 2 
or 0(1) in d > 3. 



7 Harness process in continuous time 

The harness process in continuous time can be constructed as the L — > oo 
limit of the harness process with random sequential update, defined like the 
sub-lattice parallel dynamics but with ([2]) replaced by 

p(d^i^ r - 1 ) = E c 2lVl Vl,J ■ 

3=0 

.\{S(hl - hr^UdhJ I 'norm. (74) 

i=tj i 

The time t for the Poisson clocks of rate one in the harness process in continuous 
time is related to the microscopic time r through t = t/L. The measure (|T|) is 
also invariant under the dynamics (|74l) . and we still take it as initial condition. 

Here we give numerical results for this model, indicating that the asymptotic 
forms (HU (|15[) (|16[) in Proposition [T] may still be valid, with some rescaling. 
The initial condition is drawn using the Fourier modes, which are independent 
Gaussian variables. We then run the dynamics for a time t% + 1 and measure 
the correlations: 

9^(t,j) = jjr f eW 2 - CO - ^ +e> ) (75) 

i=0 1 t'=0 

fl^fci) = |E f E^'^ - ^')(^' +t+1) - (76) 

1=1 1 t'=0 
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i=l 1 t'=Q 

The results are displayed in Fig. 2 and Fig. 3, where the upper indices L, t\ have 
been omitted for clarity, while the values of L and t\ appear in the captions. 

Fig. 2 shows relaxation as function of time, with both the numerical results 
as described above, and the corresponding exact results for the (oe: odd-even) 
sub-lattice parallel dynamics taken from Prop. [TJ The two dynamics differ for 
small time but follow similar asymptotics at large time. The function g°2(t, 1) 
is not shown on Fig. 2 because, as mentioned in Remark 1, it is proportional to 
^22 (*> 0)> an d would therefore yield the same scaled curve. 



1.6 

1.5 -K 
1.4 
1.3 
1.2 



2- 3 /2 (7 rf)l/2 5ll(M ) 

-4(2^/^3/2^0) + 



2-W 7 ViM) 

-4(70^3/2^0) 




t 

Fig. 2: Random sequential updates, scaled empirical correlation functions (?ii(t, 0), 
92%(t, 0) and gu{t, 1), average taken over space L = 10 6 and time ti = 1000, and 
scaled correlation functions g°i(t, 0), 322(^)0) of sub-lattice parallel dynamics. 



Fig. 3 shows the variation in space of the space-time correlations at a given 
large time t — 10, together with a fit inspired by l|14 p (|15 p (|16 p . 
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-0.6 1 1 1 1 1 1 

2 4 6 8 10 

3 



Fig. 3: Random sequential updates, scaled empirical correlation functions gu(t,j), 
922(t,j) and gi2 (t,j), all at time t = 10, average taken over space L = 10 6 and time 
ti — 1000, and conjectured asymptotics similar to sub-lattice parallel dynamics. 
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